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Abstract. Measurements of double helicity asymmetries for inclusive hadron production in polar- 
ized proton-proton collisions are sensitive to spin-dependent parton distribution functions, in par- 
' ticular to the gluon distribution, Ag. This study presents the double helicity asymmetry and the cross 

\ section for tj production (p + p — > 77 +X). The cross section measurement yields valuable data for 

the extraction of the fragmentation functions, which are unknown at present but are a prerequisite 
7— I ' for the extraction of the spin-dependent parton distribution functions. 

> ' 

! ' Keywords: helicity asymmetry, gluon polarization 

^ ! PACS: 13.85.Ni, 13.88.+e, 21.10.Hw, 25.40.Ep 

o : 

INTRODUCTION 

\0 • 
O ' 

; Present knowledge about spin-dependent parton distribution functions (PDFs) in the 

nucleon mainly comes from next-to-leading order (NLO) QCD fits (see, e.g., 
to the spin-dependent structure function g\ as measured in polarized inclusive deep- 
<j_) \ inelastic scattering (DIS) experiments (see, e.g., H, 0]). The resulting spin-dependent 

gluon PDFs have rather large uncertainties due to the fact that the exchanged virtual 
photon does not couple directly, i.e., at leading order, to the gluon. Thus, additional 
data from polarized pp scattering in which longitudinally polarized gluons are directly 
probed via scattering off longitudinally polarized gluons or quarks should greatly reduce 
the uncertainties in the NLO fits, especially for the gluon distribution. 

First results from the PHENIX and the STAR collaborations at the Relativistic Heavy 
Ion Collider (RHIC) are available, measuring double helicity asymmetries in inclusive 
7t° and jet production [@], respectively. The double helicity asymmetry is defined as 

A LL = a+ + ~ a+ + with Ac°cYAf a ®Af b ®Ac ab ^ cX '®D h c (1) 

where the cross section a ++ (d + ~) describes the reaction where both protons have the 
same (opposite) helicity. The spin-dependent term is given on the right hand side, where 
Af a , Afi, represent the spin-dependent PDFs for quarks (u,d,s) and gluons, and Ad 
are the spin-dependent hard scattering cross sections calculable in perturbative QCD. 
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The fragmentation functions (FFs) D h c represent the probability for a certain parton c 
to fragment into a certain hadron h, and thus they are not needed in the case of jet 
production. It is apparent that the extraction of the spin-dependent PDFs from, e.g., 
the PHENIX 7t° data, critically depends on the knowledge of the FFs. It therefore is 
desirable to extract the spin-dependent PDFs using different channels, not only because 
of possibly different systematic uncertainties involved in the All measurement, but also 
because of the different FFs involved. The latter can be extracted based on cross section 
measurements from independent data, mainly from e + e~ colliders. 

This study focuses on All m 1 production (p + p — ► r\ +X). In contrast to the n° 
the r\ FFs are unknown, but at least the quark FFs can be reasonably well constrained 
based on the existing r\ cross section measurements from e + e~ collider data. The 
extraction of gluon FFs requires either e + e~ data taken in a wide range of y/s values, 
or cross section measurements from, e.g., pp scattering. The r\ cross section has been 
recently extracted by the PHENIX collaboration based on the data taken in 2003 0]. 
An extraction based on the data from 2005, with improved statistics and a wider range 
in transverse momentum is in progress. It will serve as an important input for the r\ FF 
extraction which is in preparation |^]. 

Apart from serving as input in a global fit in order to constrain the spin-dependent 
gluon PDF, the All m f] production might have a second interesting application. Due to 
the 5-quark content in the r\ wave function which is absent in the % , it is expected that 
the fragmentation from s-quarks into 77 mesons is larger than the one into n° mesons. 
If the s-quark contribution turns out to be significant, the All in f] production might 
help to constrain the spin-dependent s-quark PDF, either in a global fit or by a dedicated 
analysis focusing on possible differences between the All in 7T° and r\ production. 

RECONSTRUCTION OF 77 MESONS 

The 77 meson is reconstructed via its main decay channel 77 — > 77 with a branching ratio 
of about 40%. The three-body decay 77 — ► 7t + n~n° not only has a smaller branching 
ratio of about 23%, but in addition also has a smaller acceptance in the PHENIX 
spectrometer and therefore has not been considered yet. The data were taken at the 
PHENIX [9] experiment in 2005. The primary detector used in this analysis is the 
electromagnetic calorimeter (EmCal), located at a radial distance of about 5 m from 
the beam pipe. It covers the pseudo-rapidity range \r\\ < 0.35 and has an azimuthal 
acceptance of = %. The EmCal consists of eight sectors, six of which are composed 
of a total of 15552 lead-scintillator (PbSc) sandwich modules (5.5 cm x 5.5 cm x 
37.5 cm), and two sectors of lead-glass (PbGl) cherenkov calorimeters, consisting of 
a total of 9216 modules (4 cm x 4 cm x 40 cm). 

A cluster in the EmCal is assumed to originate from a photon if the following criteria 
are met. First, the cluster may not be associated with a signal from a charged particle in 
the Pad Chamber just in front of the EmCal (PC3); an exception is made if the signal 
position in the EmCal and in the PC3 are aligned in such a way that the particle likely 
came from the vertex on a straight line, i.e., it was not bent in the central magnetic 
field. In this case the cluster is accepted as a photon candidate since it is assumed that 
the original photon converted into an e + e pair just before the PC3. Furthermore, since 



electromagnetic showers in the EmCal are not confined within a single module, a shower 
profile analysis can be used in order to reject hadrons which usually produce broader 
showers than photons. Since hadrons are slower than photons, an additional time of 
flight cut is used for the photon identification. 

Using all possible pairs of photon candidates, the two-photon invariant mass can be 
calculated. In order to exclude clusters with potentially wrongly reconstructed energies 
from the calculation, the module with the largest energy deposition in a cluster may 
not be in the outermost two columns or rows of an EmCal sector, and there may not 
be a noisy or dead module in the eight surrounding modules. An upper limit of 0.7 
is placed on the energy asymmetry [E\ —E2)/(E\ +E%) of the two cluster energies 
E\ and E2 m order to reduce the combinatorial background. Finally, the transverse 
momentum pj of the di-photon is required to be larger than 2 GeV/c. The invariant mass 
distribution for two different pj bins is shown in Fig[TJ In the vicinity of the r\ peak the 
distribution is well described by a fit to a Gaussian plus a second-order polynomial. For 
the asymmetry calculation the events within a 2a window around the mean of the 77 
peak are used, whereby the number of background events N BG is given by the integral 
under the polynomial within this window. 
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FIGURE 1. Two photon invariant mass distribution for two different pj bins. 



DOUBLE HELICITY ASYMMETRY 



Experimentally, the double helicity asymmetry (Eqn. Q]) translates into 

1 N++-RN+- , L ++ 

A LL = - — n — r — — , with R = ^±, (2) 

\P b \\Py\N ++ +RN + - L+_' 

where N++ (N+-) is the experimental yield for the case where the beams have the 
same (opposite) helicity. The relative luminosity R is measured by a coincident signal in 
two beam-beam counters, which have a full azimuthal coverage at a distance of about 
±1.4 m from the target. The achieved uncertainty on R is on the order of 10 -4 . The 
polarizations of the two colliding beams at RHIC are denoted by Pg and Py. The degree 
of polarization is determined from the combined information of a pC polarimeter lllOll . 
using an unpolarized ultra-thin carbon target, and from pp scattering, using a polarized 



atomic hydrogen gas-jet target UM- The preliminary average polarization value for the 
data from 2005 is 47% with an uncertainty of 20% per beam, leading to a 40% scale 
uncertainty in the preliminary All result. 



The double helicity asymmetry for r\ production is calculated as 
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where A L ' L is the asymmetry from the events in the 2o window around the mean of 
the r\ peak, and Af£ is the background asymmetry calculated in an invariant mass region 
below (300-400 MeV/c 2 ) and above (700-800 MeV/c 2 ) the tj peak (see Fig. Q). The 
asymmetries in the two regions are consistent with each other and therefore have been 
combined to form Af^. The ratio r gives the number of background events divided by 
all events in the 2o window around the mean of the 77 peak. The background corrected 
asymmetry A^ L as a function of pj is shown in Fig.[2l As discussed above, the asymme- 
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FIGURE 2. Double helicity asymmetry for the inclusive Tj production as a function of pj. 

try can only be related to the spin-dependent gluon distribution after the r\ fragmentation 
functions have been extracted. Similar to the PHENIX results on the n° asymmetry, the 
7] asymmetry is consistent with zero over the measured range. It should be noted that 
the data taken in 2006 will improve the statistical error for A\ L by about a factor of three. 
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